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Abstract. Discrete-time quantum walks (QWs) represent robust and versatile platforms for 
the controlled engineering of single particle quantum dynamics, and have attracted special 
attention due to their algorithmic applications in quantum information science. Even in their 
simplest ID architectures, they display complex topological phenomena, which can be employed 
in the systematic study of topological quantum phase transitions 1 . Due to the exponential 
scaling in the number of resources required, most experimental realizations of QWs up to date 
have been limited to single particles, with only a few implementations involving correlated 
quantum pairs. In this article we study applications of quantum walks in the controlled 
dynamical engineering of entanglement in bipartite bosonic systems. We show that quantum 
walks can be employed in the transition from mode entanglement, where indistinguishability of 
the quantum particles plays a key role, to the standard type of entanglement associated with 
distinguishable particles. We also show that, by carefully tailoring the steps in the QWs, as 
well as the initial state for the quantum walker, it is possible to preserve the entanglement 
content by topological protection. The underlying mechanism that allows for the possibility of 
both entanglement engineering and entanglement protection is the strong "spin-orbit" coupling 
induced by the QW. We anticipate that the results reported here can be employed for the 
controlled emulation of quantum correlations in topological phases. 



1. Introduction 

Discrete-time quantum walks (QWs), originally inspired by the Feynman path-integral approach 
for computing observables [2], are the quantum extension of classical random walks, where clas- 
sical probabilities are replaced by quantum imaginary amplitudes, so that the multiple paths 
of the quantum walker can interfere, a feature which results in the probability of the quantum 
walker being quite distinct from their classical counterpart [3]. Within the realm of classical 
statistical physics, random walks have been extensively employed as fundamental models de- 
scribing the complex dynamics of several systems of interest, ranging from transport properties 
in biological and material compounds, to the dynamical evolution of financial markets. On the 
other hand, the analogous quantum walks have been used to model a large number of dynamical 
physical processes involving quantum coherence, such as photosynthesis [4J, diffusion [5j, and 
vortex transport [6J. In recent years, QWs have attracted significant attention as they provide 
a universal platform for the study of quantum information protocols, allowing for speed-up in 
search algorithms 018], and for universal simulation of quantum circuits [9J. Moreover, it has 
recently become apparent that QWs can be employed for the efficient dynamical exploration of 



a wide range of non-trivial topological phases in ID and 2D jTJ [TUJ • In particular, the simplest 
ID random walk that presents non-trivial topological features is analogous to the SSH model in 
polyacetylene [TT] . Several experimental realizations of QWs have readily been accomplished, 
employing trapped ions [12], ultra-cold atoms in optical lattices [13], and photons [14J. The lat- 
ter are of particular interest, as light particles can be easily produced, manipulated and detected. 
However, due to the exponential growth in the Hilbert space size and in the number resources, 
only a few QW implementations up to date involve correlated quantum particles [15} [TT]. 

In this paper, we analyze the use of QWs in the controlled engineering of entanglement be- 
tween bosonic quantum particles. We discuss in detail the different types of entanglement that 
can be observed between the quantum particles, namely mode entanglement of indistinguish- 
able particles, and standard entanglement between separate distinct particles, the latter being a 
fundamental resource for quantum communications. The entanglement of elementary states of 
indistinguishable particles was discussed in detail in the classic book by A. Peres [16J, and has 
been the subject of continued active research [T8J [19] . Here we show that quantum walks can be 
employed to engineer controlled transitions between different types of entanglement. This can 
be explained in terms of coupling between different internal degrees of freedom in the quantum 
walker, such as spin and spatial or orbital degrees of freedom, a result which can be exploited 
for the controlled dynamical engineering of entanglement. This form of "spin-orbit" coupling 
mechanism induced by the quantum walk is also at the heart of topological phenomena, as first 
introduced in the seminal work by Kitagawa and co-workers [1] . 

Here we demonstrate that by carefully engineering the parameters in the quantum walk, such 
as the quantum coin, and the initial state of the quantum walker, it is possible to preserve the 
entanglement content from noise and perturbations, a feature which can, in turn, be explained in 
terms of topological protection. The term "topological protection" , usually refers to immunity 
to local noise and back-scattering of highly symmetrical ground state wave- functions of quantum 
systems characterized by some form of topological order. Usually, such topological bound states 
display strong spin-orbit coupling and are protected against perturbations by time-reversal sym- 
metry. The idea of entanglement preservation by topological protection, is to prepare the initial 
state of the system in such a way that there is a significant overlap with the topological bound 
states. Here, we show that by tailoring the "band-structure" characterizing the QW, this can 
be done by choosing particular values for the quantum coin, it is possible to simulate effective 
Hamiltonians which are gapless, showing a Dirac-like dispersion relation. In the vicinity of such 
Dirac points the emergence of robust localized bound states is expected. Therefore, by encoding 
the initial state of the bipartite quantum walker in the topologically protected bound state basis, 
it possible to protect its entanglement content. 

The article is structured as follows, in Section II we review the formalism for quantum walks 
of single and composite systems, as well as the relevant entanglement measures. Next in Section 
III, we present our results for entanglement engineering in systems of indistinguishable particles. 
In Section IV, we present application of quantum walks in topological protection. Finally, in 
Section V we present the conclusions. 



2. Formalism 

2.1. Basic Formalism: QW for single 2-level systems in ID 

The simplest quantum walk consists of a protocol for a 2-level particle (t,i), moving in dimension. 
The unitary evolution operator U for a single step in the quantum walk is given by: 



U = TRv(0), 



(1) 



where f = Y, x I t)(t I ® \x + l)(x\ + | \ ® \x — l)(x\, is a spin-dependent translation by 
one site x, R${9) — exp (^—i9v • <r/2^ is the coin operator given by an arbitrary SU(2) rotation 

around an unitary vector v in the Bloch sphere, and c — ((j^, a y ,a z ) are the Pauli matrices. It 
is the quantum equivalent of randomly tossing a coin and choosing which way the particle will 
move. Conventionally, the rotation is set to R y — e - l0cr y/ 2 . The unitary evolution operator U 
can be considered as a stroboscopic regularization of a static effective Hamiltonian H e ^ given 
by: 

U = e - iA * G . (2) 

We note that the spin-dependent translation can be written in quasi-momentum space as 
T = \^dke %k(Jz ® \k)(k\, which readily reveals the spin-orbit coupling mechanism which is at 
the heart of the entanglement engineering and topological protection mechanisms described in 
this paper. In quasi-momentum space, the effective Hamiltonian up to a constant energy offset, 
can be written as: 

# eff = r dk[E(k)n(k) • ff\ ® \k)(k\. (3) 

J — TV 

Now n{k) defines the quantization axis for the spinor eigenstates at each quasi-momentum 
k. For a coin operator different from the identity operator, the quasi-energy is given by 
cos(E(k)) = cos(#/2) cos(fc), typically corresponding to a 2-band structure, for a 2-level sys- 
tem. The quasi-energy spectrum is plotted in Figure 1 for different values of the coin parameter 
6. For |0| > the spectrum has a gap, and for 9 w 0, 2tv the spectrum becomes gapless, and 
presents a Dirac-like dispersion relation, where the quasi-energy vs quasi-momentum are linearly 
related. This already explains why QWs in their simplest architectures can be used to efficiently 
simulate dynamical properties of complex condensed matter system with topological properties. 




Figure 1. Band structures for the simplest ID QW in a 2-level system. By tuning the quantum 
coin parameter 9 it is possible simulate a gapped system, or a gapless system with Dirac points. 
9 = (black), 9 = tt/2 (red), 9 = tt (green), 9 = 3n/2 (orange), 9 = 2tt (blue). 



In general, non-trivial topological character requires at least two independent parameters, so 
that an adiabatic variation of these parameters leads to the accumulation of a Berry phase [20] . 



This can be done by introducing a second rotation R y (d 2 ), which is called in the literature the 
"split-step" quantum walk [I]. 

2.2. Split- step quantum walk 

The split-step QW protocol reads : 

17(01, 6 2 ) = f i R y {e 2 )f^R y {e 1 ). (4) 

where: 

T-r = El t)(t \®\x + l)(x\ + \l)(l\®\x)(x\, (5) 

X 

Ti = EU)(4-l®|a:-l)^l + lt>(tl®k)(a:|- (6) 

X 

It has been shown analytically [lj that this system still presents chiral symmetry, meaning 
that the positions {n(k)} of the eigenstates of U (61,62) lye on a circle of the Bloch sphere. 
Nevertheless, depending on the values of 6\ and 62, the set {n(k)} can describe the full circle 
(W = 1), or not (W = 0), where W is the winding number associated to the topological struc- 
ture of the spectrum. The spectrum presents a gap in both cases |tan(02/2)/ tan(0i/2)| < 1 
(W = 1) and |tan(0 2 /2)/tan(0i/2)| > 1 (W = 0). When | tan(0 2 /2)/tan(0i/2)| = 1, the gap 
closes (at quasi-energy or 7r), and one can not properly define a winding number. In this case, 
the corresponding 0(resp. tt) quasi-energy state corresponds to a so called topological bound 
state. The topological phase diagram with two coin parameters becomes a mesh-grid of alter- 
nating squares, characterized by alternating topological numbers W = 0, 1 (see fig. (2)). We 
can notice that when 62 = 0, we recover the standard quantum walk protocol operator. 




Figure 2. Diagram representing the topological nature of the split-step QW. Full red (resp. 
dashed black) lines represent configurations that give a gapless spectrum where the gap closes 
at quasi-energy (resp. tt). The numbers (0 and 1) written in the different zone are the winding 
numbers of the given choice of angles (0i, 62). 



In order to experimentally observe the presence of the topological nature of the split-step 
QW, one can decide to make one of the rotations inhomogeneous; for instance, make 62 vary as 
a function of x, going from 0^" (when x —00) to 0^ (when x — > +00). If one line where the 



gap closes is crossed by the trajectory of (#1, #2), we know that the local spectra are gapped 
far from the crossing point, and gapless in the vicinity of this point (see fig. (3)). The as- 
sociated bound state (say of quasi-energy) will be localized around the zone where the (#1, 
62) trajectory crosses the line (say x = 0). Already at a few sites from x — 0, there is no 
arbitrary small quasi-energy states that could permit the propagation of the quasi-energy 
component of the wavepacket, and consequently the overlap between the initial wavepacket and 
the localized bound state is topologically "trapped" around x — 0. Then, after propagation, the 
probability to detect the particle localized around x — remains significant compared to the 
case in which bound states are not present in the spectrum. Indeed, the basic feature of the 
QW evolution of an initially localized (around x — 0) particle is a ballistic expansion which is 
in general anisotropic (not symmetric in ID) and where the probability for the particle to stay 
at its initial position decreases very quickly. This difference of behaviors has been indeed exper- 
imentally observed [23], and clearly shows the important impact that the topological structure 
of the system has on the dynamics. Note that the trajectory (#1, 62) has to keep 9\ uniform 
in order to preserve the chiral symmetry of the system, but the variations of (#2) can be arbitrary. 
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Figure 3. Spatially varying coin parameter #2(2?) characterizing the inhomogenous split-step 
quantum walk. The horizontal axis represent the spatial degrees of freedom for the ID QW (re- 
direction). This inhomogeneous step allows to assign different values to the parameter across 
topological sectors characterized by different winding numbers (W = 0, 1) and in this way create 
a boundary. At the spatial location of the topological boundary (x = 0) a topologically protected 
bound state is expected. This feature can be used for entanglement topological protection, by 
encoding the initial state of the quantum walker in a localized state in the vicinity of this 
boundary. 



2.3. Basic Formalism: QW for two identical 2-level systems in ID 

We consider now the dynamics of two identical bosonic particles, such as linearly polarized 
photons, propagating in the same QW protocol, represented by identical unitary evolution op- 
erators U(0ii #2), as defined in the previous subsection. Since photons do not interact, we focus 
on the evolution of entanglement on the two-particle dynamics in order to analyze to impact of 
quantum correlations in the system. 

The Hilbert space of a bipartite system is given by H = Ha ® Wb, and the step operator 
is now U = Ua ® Ub- In the case of indistinguishable particles, as considered in this work, 
the structure of the Hilbert space is different because both particles "live" in the same Fock 



space. The adapted formalism in this context is the second quantization, that automatically 
symmetrizes (resp. antisymmetrizes) the bosonic (resp. fermionic) wavefunctions. 

Let us first consider the case that the two photons are at the same spatial location x, can 
be described by identical frequencies (i.e., monochromatic case), and can only be labeled by 
the polarization degree of freedom (modes f and \). Any bosonic two-particle wavefunction 
can be built by the superposition of the states of the canonical symmetric basis composed by 
a^a^\mathrmvac} = | tt) > ajajjvac) = | H) and ajajjvac) = ^ (| It) + I ti)) where |vac) 
represents the vacuum state of the Fock space. These states can be written in a more syn- 
thetic form, in terms of mode (t or j.) with an occupation number : a|a||vac) = \2 *f 5 \), 

ajajjvac) = |0 t?2 ^), ajajjvac) = |1 "f 5 1 4-)- Additionally, another advantage of this notation 
is that it clearly appears now that the entanglement in polarization of ^ (| J/f) + | f^)) is just 
the consequence of the bosonic character of the photon pair because in the second quantization 
language, it is written as a product state in modes |1 ^ , 1 ^), as well as the two other states 
mentioned above. In particular, the state ^ (| ft) + I ID) which is maximally entangled in 

the polarization representation can be rewritten (|2 f, \) + |0 t, 2 \)) in the mode repre- 
sentation, which shows that this state is maximally entangled in both mode and polarization 
representations fT8| [T9]. 

In order to respect the global symmetry of exchange of bosonic particles, the states describ- 
ing the case in which one particle is at position x\ and the other at X2 are: on one hand, any 
possible tensorial product of the symmetric state in positions (4= (\xi,X2) + \x2,xi))) and one 

of the symmetric state in spins (| ft}? I 11) > or ^ (I it) + I tt))) or > on the other hand, the 
tensorial product of the antisymmetric state in positions (\xi,X2) — \x2,xi))) and the anti- 
symmetric state in spins (4^ (| ^t) — I tt)))- Since each photon evolves with the same unitary 
evolution operator, an initially symmetric state with respect to the exchange of particles will 
remain symmetric after an arbitrary evolution time. We note that in the case of photon pairs 
produced by Type-II Spontaneous Parametric Down Conversion (SPDC), the two-photon state 
has the general form = J dq\ J dq 2 J dw\ J dw2[(/>Hv(gi^i]Q2^2)\Hi,qi,wi]V2,q2,W2) + 
^vh(qi^i] Q2, ^2)|Vi, qi : wi] H2, q2, ^2)]? where (i7, V) stand for orthogonal linear polarizations 
(t,t), and (q,w) represent the spatial and frequency degrees of freedom, respectively. For the 
case of indistinguishable particles, (i.e., no temporal or spatial walk off), the two-photon ampli- 
tude wave-functions should be identical 4>vh(qi, wi] ^2) = feyfc^i;^)^)^ which forces 
the spatial and frequency part of the wave-function to be symmetric. In this case, the state is 
maximally entangled in polarization (i.e., Bell state). Any uncontrolled spatial or temporal walk 
off reduces the indistinguishability of the photon-pairs, and therefore the entanglement content. 

In the following, in order to quantify the entanglement of a given two-particle state, we 
will compute the logarithmic negativity either in mode representation, when both particles are 
localized at the same position, or in polarization representation when the particles are physically 
separated (i.e., the overlap between the two wavepackets can be neglected), and can then be 
treated as distinguishable particles. 

3. Entanglement engineering by discrete-time quantum walks 

In order to illustrate the transversality of the concept of entanglement between identical particles, 
we engineer a split-step QW protocol for which two indistinguishable photons initially localized 
around x = can be measured with a reasonably high probability at different positions after 



the propagation, and can then be distinguished. With a well choosen set of coin parameters 
6^), we will show that an initially mode-entangled state leads to non- vanishing entanglement 
in polarization between the two photons, whereas an initial state which is separable in terms 
of modes does not exhibit this behavior. We will refer to this phenomenon by "entanglement 
conversion" , from mode entanglement to polarization entanglement, during a temporal transition 
from indistinguishable to distinguishable particles. 



3.1. Numerical Results 

We numerically compute the split-step QW dynamics of the two following initial states: 
|* A > = \xi = 0,x 2 = 0) ® ^ (I It) + I tl» and \V B ) = \x x = 0,x 2 = 0) ® ^ (| ft) + I II))- 

The choosen coin parameters are 9\ = | and 9f = ±|, and then the topological aspects of the 
split-step QW are not involved here since the whole system belongs to the same topological zone 
(see fig. (2)). 
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Figure 4. Spatial probability distributions in logarithmic scale (color online) of \^a) (top) and 
\^b) (bottom) after 60 iterations. 



We clearly see on figure (4) that there is almost no difference on the probability distributions 
between the two cases. Indeed, both present important symmetric peaks of probability on the 
off-diagonal line : x\ + x 2 = 0, which represents one photon on each side. For each initial state 
and at each time step n, we first project the two-particle state (lJ(9i,9f) J7(0i, 9f)^j \^a,b) 
on the two-position state |x max , — x max ) corresponding to the maximum probability lying on 
the off-diagonal line, then trace out the position degrees of freedom, and finally normalize the 
result. This procedure gives a two-polarizations vector |^o), called "projected" state. Hence, 
one can build the 4x4 "reduced" density matrix p = |V ; o)(V ; o| an d evaluate the entanglement 
between polarizations by computing the polarization Negativity: N(p) = Tr(|p|) — 1, where p is 
the partially transposed density operator. 

A radical difference is observed in terms of polarization entanglement during time evolution. 
In the case of a non mode-entangled initial state (|^^)), some fluctuations of the polarization 
negativity are observed, close to the minimum value 0. On the other hand, if the state is 
initially entangled in modes (|\I/#)), part of this entanglement is transferred (or converted) to a 
traditional polarization entanglement during time evolution. 

3.2. Discussion 

The split-step QW protocol allowed us to propose a simple realization exhibiting a dynamical 
transition from indistinguishable to distinguishable bipartite system. Additionally, it can offer 
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Figure 5. The red solid (resp. dashed blue) line represents the computed polarization 
Negativity of the measured reduced density matrix for initial state \^a) (resp. I^b)), versus 
time n. In both cases, it is maximum (N(p) = 1) at the initial time. 



a nearly direct observation of mode entanglement between identical particles. Let us insist on 
the fact that at short {i.e. experimentally meaningful) times, the measurability of this effect 
is enhanced because the probability to measure the two photons on the anti-diagonal line is 
decreasing with time. 

In the present section, we have avoided the presence of a non-trivial topological structure of 
the split-step QW in order to focus on the nature of entanglement between identical particles. 
We will now include this important feature, and focus on the interplay between entanglement 
and the eventual presence of topological bound states in the spectrum of the single-particle 
evolution operator. 



4. Entanglement protection by discrete-time quantum walks 

Let us consider a configuration of Of) that crosses one (fig. (6), left panel) or two (fig. (6), 
right panel) gap closing line(s). 

Here, we show that when an entangled two-particle state (\^b)) follows a QW protocol for 
which there is one bound state |^bs) localized in x — 0, the two-particle wavepacket presents 
an important localized part around the origin x\ — X2 — 0, which can be locally identified to 
I^bs) ® I^bs) (see fig. (7), top panel). When two gap closing lines are crossed, the localization 
phenomenon around the origin is even stronger (see fig. (7), bottom panel). In this case, two 
different bound states exist in the spectrum of the single particle evolution operator, one l^s) 
of quasi-energy, and the other |^§g) of tv quasi-energy. 

We suppose that both photons are now detected at the origin, and compare the mode 
entanglement during time evolution. This mathematically corresponds to projecting the two 
particle state on \x± = 0,X2 = 0), building up the density matrix represented in modes p m , and 
compute the mode Negativity N(p m ). 




Figure 6. Configurations of the coin parameters leading to the presence of one (left) or two 
(right) topological bound states. The red lines describe the trajectories (#1, #2 (a?), x G 5R). 
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Figure 7. Spatial probability distributions in logarithmic scale (color online) of \^b) after 60 
iterations. The spectrum of the single-particle evolution operator contains one single (top) or 
two (bottom) bound states. 



We can easily observe that when two bound states are present in the system, the phenomenon 
of localization around the origin is clearly enhanced, and thus the probability of detecting both 
particles in x — is increased. We can also notice that the entanglement vanishes quickly 
during time evolution in the first case (1 bound state present), whereas in the other scenario (2 
bound states present), the entanglement oscillates between 0.2 and 0.4 without decreasing at all. 



4.1. Discussion 

Similarly to the single particle case, the main feature of the QW is a ballistic expansion of an 
initially localized wavepacket. Here, after a certain time, an important part of the wavefunction 
leaves the origin via ballistic expansion, but what remains in the center is the overlap of the 
initial state with the bound state(s). In the first case, we can approximate this fraction of 
the state by a product state: ~ Xl ,x 2 ttO |*bs) ® |^bs) 5 where |^bs) represents the zero 
energy bound state present in the first case. We can then easily understand the fact that the 
entanglement rapidly decreases because of the presence of one single bound state. 

In the second column case, two bound states are present (|^bs) an( ^ I^bs))> an( ^ as previously 
explained, after a certain time of propagation, the central part of the wavefunction is mainly 
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Figure 8. Computed Negativity of mode entanglement when the two photons are detected at 
the origin, for the split-step QW evolution of the initial state \^b)- The red full line corresponds 
to the presence of a single bound state, and the oscillatory behavior displayed in the blue dashed 
line is a clear signature of to the presence of two bound states with quasi-energy difference tt. 



due to the presence of bound states. When two bound states (or more) are present, we can not 
conclude that the form of the state is a product: in general, it would be a coherent superposi- 
tion of the different bound states. This qualitatively explains the important difference observed 
in the evolutions of entanglement. In the latter case, the non vanishing negativity is directly 
related to the topological structure of the system and thus, this entanglement should be robust 
to the presence of impurities, or disorder. 

The co-existence of pairs of bound states with quasi-energy difference tt is a robust 
phenomenon with topological origin, characteristic of periodically driven systems |21j . In our 
case we find a similar effect for a static effective Hamiltonian due to the presence of two distinct 
bound states with quasi-energy E — and E — 7r, when crossing two distinct topological 
boundaries for our specific choice of coin operators. The co-existence of these bound states, 
in turn, explains the periodic oscillations in the Logarithmic Negativity displayed in Figure 8, 
and represents a novel feature that can be exploited for entanglement protection. The insight 
into topological protection gained by this approach can help to explore topological properties in 
novel correlated materials, where current technology readily enables experimental realizations. 

5. Outlook 

To conclude, in this article we proposed the use quantum walks (QWs) for the controlled engi- 
neering and protection of entanglement in bosonic bipartite systems. The underlying mechanism 
that allows for manipulation and control of entanglement, as well as for its topological protec- 
tion, is the strong coupling between internal degrees of freedom in the quantum walker, which 
plays an analogous role to "spin-orbit" coupling in topological systems. In particular, we showed 
that QWs can be employed in the transition from bound entanglement to mode entanglement 
and that such entanglement content can be topologically protected by tailoring the quantum 
coin, and the initial state for the quantum walker. The extension of QWs to correlated particles 
pairs opens the door to simulations of quantum correlations in condensed matter systems with 



topological order. 



The results presented here have direct experimental applications in photonic quantum 
systems, and could improve the performance of photonic quantum technologies [22J. In this 
regard, we note that the phase stability of multi-path QW interferometers is not expected to 
represent a significant experimental limitation, since topologically protected states are expected 
to be resilient against local phase noise and perturbations. However, temporal or spatial walk- 
off in the system with increased number of iterations can play a deleterious dephasing effect 
which, in combination with losses, can limit the number of iterations. Current bulk optics 
implementations are limited to n — 7 iterations at the most [23] . 
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